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Abstract 


We find a solution of an unknown time-dependent diffusivity a(t) in a lin- 
ear inverse parabolic problem by a modified genetic algorithm. At first, 
it is shown that under certain conditions of data, there exists at least 
one solution for unknown a(t) in (a(t), T(z, t)), which is a solution to the 
corresponding problem. Then, an optimal estimation for unknown a(t) is 
found by applying the least-squares method and a modified genetic algo- 
rithm. Results show that an excellent estimation can be obtained by the 
implementation of a modified real-valued genetic algorithm within an Intel 
Pentium (R) dual-core CPU with a clock speed of 2.4 GHz. 
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1 Introduction 


This article considers helium diffusion kinetics, which is important for materi- 
als in which helium measurements are made, particularly for thermochronol- 
ogy (see [2, 4, 3]). Solutions of an inverse problem for the helium production- 
diffusion equation are explored in this research article. This solution is a 
time-temperature path derived from apatite fission-track length distribution 
to characterize the response of apatite helium ages to thermal histories in- 
volving partial retention. Mathematically, the helium production-diffusion 
equation without any source of radiogenic production in a rectangle domain 
Q = {(x,y):0<a4<1,0<y<1} can be written as (see [2, 4, 3]) 


OT OT C&T 
a — (Naa + az (z,y) € QO, t > 0, (la) 
T(x, y,0) = f(a,y), (x,y) € Q, (1b) 
T(0,y,t) = go(y,t), O<ys lt 20, (1c) 
TU, = nly, t), O<y<1,t2>0, (1d) 
oF (0, 0,t = ho(z,t), O0<a<1,t>0, (le) 
y 
oT (edt = hi(a,t), O0<a<1,t>0, (1f) 
y 


where f(z, y), go(y,t), g1(y,t), ho(a,t), and hi(#,t) are known functions and 
a(t) is the time-dependent diffusivity and should be determined in this inverse 
problem. 

Here the special case ho = h, = 0 is considered, that is, 


OT OT &T 
at _ a(t)( Ox? y Oy? iF (x,y) € Q, t> 0, (2a) 
T(x,y,0) = f(x,y), (x,y) €Q, (2b) 
T(0,y,t) = go(y, t), O<y< 1,t>0, (2c) 
T(1,y,t) = gly, t), O<y<1,t20, (2d) 
OT OT 
— t)= — 1,t)= Saree Let 
aye ) By ,t) =0, 0<2r<1,t>0, (2e) 


together with an overspecified condition for fixed points 29, yo € (0,1), as 
Cannon considered in [7, p. 191], 


a(t) <—(0o.vort) = p(t) (#2 0). (3) 


Problem (1) in one-dimensional case arises from the heat conduction in a ma- 
terial that is eroded or damaged. This erosion can be caused, for example, 
by radioactive. Duo to the erosion, thermal characteristics of the material 
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including its heat capacity, diffusivity, and conductivity change based on 
the level of the erosion. These characteristics can be related to time, where 
a(t) is the thermal diffusivity, which is the ratio between the thermal con- 
ductivity and the heat capacity. Lesnic [14] investigated an inverse problem 
for determining the temperature and the time-dependent thermal diffusivity 
from various additional nonlocal information (see also [7, Chapter 13]). It 
should be remarked that the inverse problem (1) with a(t) = 1 has been 
investigated by several authors; see, for example, [5, 20, 8, 22, 23, 9] and for 
a two-dimensional case [21]. Since solving inverse diffusivity problems using 
evolutionary algorithms has not been investigated yet, the solution to these 
problems using a modified genetic algorithm is considered in this article in 
the hope that this effort triggers other authors’ motivation to examine new 
artificial intelligence techniques. 


2 General solution in term of Green’s functions 


For 2 = (0,1) x (0,1) and 2 = [0,1] x [0,1], take 


5; ={x=(2,y): 2 € {0,1},0<y< 1}, 
Sy ={x=(r,y):y€ {0,1},0<a< 1}. 


Then 02 = SU Sg. Consider the two-dimensional heat conduction problem 
with nonhomogeneous boundary conditions 


pp t) =VPT (xt), xeEO,t> 0, (4a) 
T(x, 0) = f(x), x € 0, (4b) 
T (x,t) = 9(x,t), x€ $1, t>0, (4c) 
5,08) = h(x,t), x € So,t>0, (4d) 


Here go, gi, ho, hi, and f are known functions. It is known that (4) possesses 
a unique solution provided the given initial and boundary data are continuous 
functions; see, for example, [7, 10, 19]. The aim of this section is to find 
a representation of the solution of (4). To this end, as in [7], define, for 
xeER,t>0, 
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K(a,t) = exp(——). O(a, t) = = K («+ 2m,t), 


m=—Cco 


1 
VAnt 
and define the function ®(x,t) = K(a,t)K(y,t) for x = (a, y) € R?, t > 0 by 


1 exp( Ir) 
Ant 4t / 


It is well known from [12, Chapter 7], that 


®(x,t) = 


T(x,t)=Oax f= | Ox-Et)f(Edé (xeER’,t>0),  ( 
R2 


OU 
nN 


where T(x, ¢) is a solution for a pure initial value problem in the equation 


OT 9 
op & t) = V°T (x,t), 
T(x, 0) = f(x), 


for x € R? and t > 0, provided that the function f(x) is bounded and also 
(piecewise) continuous in R?. Now, we extend T(x, t) continuously to Q, 
and consider the following heat conduction problem in Q with homogeneous 
boundary conditions: 


ap Oh) = VT (et), xe 0.4 >, (6a) 
T(x, 0) = f(x), xe Q, (6b) 
T(x,t) =0, x€ 51,t>0, (6c) 

oT 

a = F x€S,,t>0. (6d) 


To solve problem (6), one can first extend f(x) to [—1,1] x [—1,1] such that 
f is odd in x and even in y. Then we extend f periodically to R? by 


f(a +2m,y+2n) = f(x,y), (n,m €Z). 


If we define, for 7,€ € R, t > 0, 


o~ (x, €,t) = A(x a & t) as O(a 7 E,t), 


then we can obtain from (5) for T(x, t) that 
Text) = [O° (eNO wn OLE naan (7 


Now, define G by 
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G(a,y, th) =O (yyy, OO" (@,6)0). (8) 


Since 6(x,t) is an even function with respect to the variable x for all values 
of t > 0, it is easy to observe for every €,7,t that 


G(x, y,t;€,n) = 0, (x,y) € Si, 


a 
79 (ts ts 80) = 0, (x,y) € So. 
Y 


So T(x, t), defined by (7), is, in fact, the solution of (6). This means that 
the function G, defined by (8), is the Green function for problem (4), so that 
the solution of (4) can be written (see [19, Chapter 6] or [13, Chapter 7]) by 


T (x,t) = [ G(x, t:€) f(@dk 
+ fae [Seat —ni8)al& en 0) 
=i ar f (mtr @)nlg Ya, (10) 


where x = (x,y), € = (€,7), € € Si in (9), and € € Sy in (10); and (0/0n;), 
in (9), denotes the directional derivative on S;. Therefore, one can write the 
solution T(a,y,t) of (4) in the following form: 


T(x, y,t) = [ S(e.u.t§. mt mdedy 
—2 I go(n, r) 2a, t—7)@*(y,n,t — r)dndr 
+ 2 | nln.) on(a —1,t—r)O*(y,n,t — 7)dndr 
eo 2f ho(n, T)0(y,t — T)O™ (a, €,t — r)d&dr 


-2 f hi(y,7T)0(y —1,t—T)O (a, €,t — r)dédr. 
Q 


Assuming that a(t) is given, one should note that under some suitable 
hypotheses on the initial and boundary conditions, the existence of a unique 
solution to (6) can be proved. 


Theorem 1. Let + > 0, 2 = (0,1) and f € L?(Q). Suppose also that 
a € L©((0,7)) and a(t) > 0, for almost every t € (0,7). Then there exists a 
unique function u of (6) satisfying 


T € L? ((0,7); Hy(Q)) NC ([0, 7]; L7(Q)) 
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and ar 
EE (0.7): Ho), 


where H}(Q) and Hy '(Q) are the usual Sobolev spaces. 


Proof. The proof of Theorem 1 follows from [16, Theorem 7.1] (see also [1, 
15, 17, 18)). 


Remark 1. If a is given and the assumptions of Theorem 1 hold, then one 
can obtain a unique solution T = T(x, y, t; a) of (6). Hence some consistency 
conditions, in boundary, should hold between a and u in (3). Consequently, 
the coefficient a is uniquely determined. 


Remark 2. One should note that Theorem 1 improves [3, Lemma 2.1]. 


3 Utilizing an appropriate transformation 


Turn back to the main problem (2), in which u(z, y, t) and a(t) are unknown. 
Denote by ¥ the family of all continuous functions a : R — R such that 
a(t) > 0, for allt € R. For a € Ff, let (see [7]) 


Then a(0) = 0 and a’(t) = a(t) > 0, and thus s = a(t) is a nonnegative 
invertible function with inverse t = 3(s). By plugging 
v(a,y, 8) = ula, y, B(s)), 
do(y, 8) = go(y, A(s)), oily, s) = my, B(s)), 
b(s) =a(B(s)), a(s) = p(8(s)), 


problem (2) is rephrased in the following form: 


Us = Vax + Vyy; (x,y) €Q, 5 >0, (11a) 

v(x, y,0) = f(a,y), (x,y) €Q, (11b) 

v(0,y, 8) = go(y, §), Uys Le So, (11c) 

v(1,y, 8) = oily, 8), O<y<1,s>0, (11d) 

vy(x,0, 8) = vy(x,1,s) =0, O0<a<1,s>0. (11le) 


In this setting, the over-specified condition (3) implies that 
b(s)vz(0, 40,8) =9(s)  (s > 0). 


y [7, Chapter 13], the solution of (11) is represented as follows: 


The time-dependent diffusion equation: An inverse diffusivity problem 39 


v(x, y,8 )= [Hem G(a,y,s;,n)d&dn 
-2 ff Mees 9eryns—Ooolncdne 02) 
+2 ff Mie 1,5— er uns— Corin, Cian. 
Applying the change of the variable ¢ = a(r) in (12), one obtains 
u(x, yt )= [| seme (x, y, a(t); €,)d&dy 
=2 ff 2 e.a(t) —atsyo*ty matt) ~ a(7)iaoins7)analr)dr 


+f aff gn, SC — 1,a(t) — a(r))O*(y, 7, a(t) — at 


For the sake of simplicity, from now on, let 


fo= fra Qd¢ = a(t) = a(r). 


Using (13) and calculating u,(x,y,t), one writes uz(x,y,t) = Ji + Jo + Js, 
where 


de [ HEMLGl~, y.alt)sé.m)dédn, 


na-af foe @, [a0 (yn, i a)go(n, 7 )dna(r)dr, 
n=2 ff Fo e-1, f aorvn f aoinrana(nar. 


It is known that 0;(a,t) = 022(x,t) (see [7]), so 


Fhe fa = we fa) = — a Soe, fo 


Since go(y,0) = gi(y,0) = 0 and @(x,0) = 0 for x # 0 (see [7]), integration 
by parts yields that 
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tpl ra) t t 
a | pote f vet tun f ago(n.r)dndr 


a8 [oe i 2 [orn [ agotn.r)| andr 


T 


tpl t ra) t 
Js =2 J O(x = 1, [ a) a forum, f a)gi(n, r)| dndr. 
T T T 


Now by applying the over-specified condition (3), we have 


eects 


7 af fo wo, fog a) |e [orn f a)90(n,7)| dndr (14) 
+2 [ [oa -uf a) [own f a)g(n,7)| andr. 


Thus one obtains the following result. 


Theorem 2. Suppose that the function f(z, y) is piecewise continuous on 2, 
that go(y,t), gi(y,t), (0/Ot)go(y,t), and (0/0t)gi(y, t) are continuous func- 
tions for 0 < y < 1, t > 0, and that go(y,0) = gi(y,0) = 0 for all y [or just 
y = yo]. Then problem (2) possesses a unique solution (u,a) if and only if 
the integral equation (14) has a unique solution a € F. 


4 Uniqueness and existence of solution for a certain 
problem 


Now, consider the problem 


aT eT PT 
a ue = (x,y) €0,t>0, (15a) 
T(x,y,0) = (z,y) €Q, (15b) 
T(0,y,t) = O<y<1,t>0, (15c) 
T(1,y,t) = vi t), O<y<1,¢20, (15d) 
ate fish Fle ht) = 0, 0<2<1,t>0, (15e) 


and the overspecified condition 


a(t) (0.0.1) =p), (#20). 


The solution of (15) is as follows: 
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T(z,y,t) = +2 [f « 9(N,T -1 2-1, f aortyn [ a)dna(r)dr. 


The integral equation (14) becomes 


a =f fa fa? a0 “on f a)g(n,7)| andr 


ff a, [a )[-a( antntan 
0/0 © 
Moreover, 


Pal —POEO) ote [an 
0 


a(t)? 


(16) 


«| er | 0) [-a(7)¥9(n,7) + 0* 22 (y, 7) 


+a(-1, f a[-atryek a(n, 1) +07 2(n,7)] dn. 


Remark 3. Since 0*(0,7,-) = 20(7,-) and 07° (0,7,-) = 26:(7,-), equation 
(16) can be rephrased in the following forms: 


y= { fa fof [@9, — a(r \o.glanar 
at a fof (09, +a (Warglndr 


Lemma 1. If g = g(t), then one can write (16) in the following form: 


7 POF Poa, f aateyar) (17) 


Proof. If g = g(t), then 


=-(0 0, t) a4 f o- 1. fof {eo +a(r)a(r) [a0(a, [ «) ~02(0, f a)]}ar 
= fa f ogee. 


Corollary 1. Under the assumption that g = g(t) is continuously differen- 
tiable and that g(0) = 0, problem (15) possesses a unique solution (wu, a) if 
and only if the nonlinear integral equation (17) possesses a unique positive 
solution a(t) that is continuous for 0 < t < T. 


i 
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Definition 1. Suppose that g = g(t) is continuously differentiable and that 
g(0) = 0. For a positive continuous function a(t), that is, a € F, define 


(La)(t) = rot aa, f ag (ear) 


Therefore, a is a solution of (17) if and only if a is a fixed point of the 
operator £, that is, £a =a. In the rest of this section, it is proved that under 
certain conditions, £ possesses a fixed point. 


For a function h defined on [0, T], let 
|h||_ = sup{|A(s)|:0<s < t}. 


For a positive constant 6, if 6 < a(t) and 6 < b(t), then, using the mean value 
theorem, one obtains 


\cay(t) — (c0y(0| << OM Fa, 6(¢—r)g'(ryar| * 


x if 6.(1,7) if a(s) — W(s)as} f(r) 


< PO (| 9 o(t — ‘(r)d a 0 —b 
5) : ¢(1,d(t-—7))g'(r)dr} |Aellella — Olletlg(é)I, 
where i a<yv< ihe b. 


Corollary 2. Suppose for some r € (0,1) that 


PsN) fac, ae- T))g' (r)dr ~6alle SP (t € (0,7). (18) 
2 0 


Then the mapping a+> £a is a contraction. 


The main result of this section will be obtained by using (17) and Corol- 
laries 1 and 2. 


Theorem 3. Suppose that p(t) and g = g(t) are continuously differentiable, 
that g(0) = 0, and that, for some r € (0,1), (18) holds. Then (15) possesses 
a unique solution (T,a). 


5 Discretization of the problem 


In this section, the following linear inverse parabolic problem is considered: 
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OT aT oF 

Ai (x,t) = a(t)( a2 7 Oye 3 Det <10<ta iw (19a) 
T(z,y,0)=f(z,y), O<z,y<1, (19b) 

T(0,y, t) = go(y, t), O<y<10<t<tm, (19c) 
T(1,y,t) = gly, ¢), O<y<1,0<t<tm, (19d) 
T(x,0,t) = ho(z, t), O0<a<10<t<ty, (19e) 
T(z,1,t)=hi(z,t), O<2<1,0<t<ty, (19f) 

with the overspecified conditions for the fixed points xo,yo € (0,1), 

T(xo,yo,t)=s(t), O<t<tm, (19g) 


where f(x,y), go(y,t), g1(y,t), ho(x, t), hi(a,t), and ¢(t) are known functions 
while a(t) and T(a,y,t) are unknown in this inverse problem. 

Now we can solve problem (19) in the least-square method. Also, we can 
consider the amount of difference between calculated values of the problem at 
the location of sensor and measured temperatures, which come from sensor as 
the cost function. When the cost function approaches to zero, the estimated 
value for a(t) approaches to exact a(t). Let 


G)= iG ay) (20) 


where T;, 7 = 1,2,3,...,m, are calculated by solving the direct heat problem. 
To do this, consider a prior guess for a(t). Also s; = s(t;), 7 =1,2,3,...,m, 
are measured temperatures. Here, when (20) reaches its minimum, the esti- 
mated value for a(t) is optimal. 


Remark 4. In this article, we use an explicit finite difference approximation 
to discretize problem (19) as follows: 


Ti j,s+1 = kas (1 Ti-1,j,s a 271 Ti, j, s { ry Tt { 1,3, s- rT; j-1,s = 2reTi, j, s 
+ reTi, 7 ae 1,s) + Lia 


FAN 1, fing Vy SS a A, 
Ti49 =flihgk), i=1,...,N—1, f=1,...,N-1, 
To,j,s =go(jk,sw), jg=l,...,.N-1, s=0,...,N—-1, 
Ti js =9(jk,sw), jg=l,...,.N—-1, s=0,...,N—-1, 
Ti,0,s =ho(ih,sw), t+=1,...,N—-—1, s=0,...,N—-1, 
T1j,s =hi(ih,sw), i=1,...,.N—1, s=0,...,N—-1, 


where x = th, y = jk, and t = sw. 
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6 Genetic algorithm 


The genetic algorithm was introduced by Holland [11] and has been used 
to solve a wide range of problems. This algorithm is a searching method 
that has been inspired by Darwinian principles of biological evolution. This 
algorithm utilizes an initial stochastic set of candidate solutions as an initial 
population. Each individual of this set is a numeric vector that estimates a 
possible solution to the problem. Then, genetic operators such as Mutation 
and Recombination are successively applied to initial populations in order to 
simulate a biological evolution. These operators push all populations toward 
optimal solutions to the problems. After a predefined number of iterations 
or termination conditions, the best individual in the population is considered 
as an optimal solution to the problem. 

In the classic genetic algorithms, the vectors are a binary string and 
present a possible solution to the problem. Indeed there are many prob- 
lems in which a real-valued vector is needed as an optimal solution. In this 
article, a real-valued genetic algorithm (RVGA) is used to estimate unknown 
a(t) in its corresponding time. The steps of our RVGA in this article are as 
follows: 


Step 1. A random set of real-valued vectors are generated as an initial pop- 
ulation. 


Step 2. Each individual in the population is evaluated by the fitness function. 


Step 3. Some individuals are selected based on their fitness by Tournament 
Selection. 


Step 4. A recombination operator is applied to the selected individuals to 
generate offsprings. 


Step 5. The Mutation operation is applied to offspring. 
Step 6. The offsprings are evaluated accordingly. 
Step 7. The population is updated. 


Step 8. Steps 3-7 are repeated, until the termination condition satisfies. 


7 A modified real-valued genetic algorithm to estimate 
unknown a(t) 


In this research article, a modified genetic algorithm is presented to determine 
unknown a(t). This algorithm employs a real-valued vector as a candidate 
solution to the problem. In this algorithm, a new step is added to the proce- 
dure of the original genetic algorithm after the Mutation operator to increase 
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the capability of the algorithm for solving the problems of this article. In 
the presented algorithm, the jth element of each vector is a real number that 
estimates a(j * .001) for 7 = 0,1,2,...,N —1. Each vector is considered as 
a chromosome in the population, and accordingly each element of chromo- 
somes is considered as a gene. Therefore, gp; is the jth gene of chromosome 
of p. To reach an acceptable estimation for a(t), equation (20) should reach 
its minimum. To meet this purpose, equation (20) is considered as the fit- 
ness function of the algorithm. When the algorithm stops, a vector with the 
lowest value of fitness is considered as the best solution of a(t). Then this 
vector is interpolated to estimate a(t) as a polynomial function. 


The steps of our presented real-valued genetic algorithm (RVGA) are as 
follows: 


Step 1. A random set of real-valued vectors is generated as the initial popu- 
lation. 


Step 2. Each individual in the population is evaluated by the fitness function. 


Step 3. Some individuals are selected based on their fitness by Tournament 
Selection. 


Step 4. To apply the recombination operator on the pair of parents, we act 
as follows: 


Jchi,j =ax Jp ,j a (1 — a) x Jp2,7> J = 1,2,3,...,N — 1, 


Jch2,j = 8X 9,; + (1-8) X Jp2,5> j= 1,2,3,...,N—1, 


Step 5. To apply the Mutation operator on the offsprings, a gene of each 
offspring is selected randomly and a small random value is added to 
it. 


Step 6. Fitness of offsprings is evaluated. 


Step 7. All other offsprings are moved toward the best offspring. Then a 
small random value is added to a gene of each offspring selected 
randomly to ensure the variety of the population. 


Step 8. The population is updated. 


Step 9. Steps 3-7 are repeated, until the termination condition satisfies. 
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Table 1: Parameters of the real-valued genetic algorithm 


Representation Real-valued vector 
Recombination Arithmetic crossover 
Recombination probability 100% 
Mutation Add a random value to one gene 
Mutation probability 1/n for each gene 
Parents selection Best 4 out of random 10 
Survivor selection All of the old generations 
replaced with new offsprings 
Number of offsprings 4 
Population size 20 
Initialization Random 
Termination condition Number of generations 


Table 1 shows the parameters of the RVGA and their roles. 


8 Numerical results and discussion 


The main purpose of this section is to show the capability of our presented 
algorithm in the previous section for solving IPP. The capability of the algo- 
rithm is investigated by an example of IPP as follows. 


Example 1. In this example, for (x,y,t) € (0,1) x (0,1) x (0,1), consider 
the IPP 


OT OT PT 


ap (tryst = aia + aa (22a) 
T(x, y,0) = (cos(my) + sin(72)), (22b) 
T(0,y,t) = ere cos(7y), (22c) 
T(1,y,t) =e-* * (cos(my) + sin(m)), (22d) 
T(x,0,t) =e"" © (sin(rx) + 1), (22e) 
T(x, 1,t) =e-" * (cos(m) + sin(m2)), (22f) 


and with the following overspecified condition 


s(tj) =7(0.1,0.1,t;)+oR, t; =(0.001)j, j =0,1,2,...,999. 
(22g) 


In this problem, the values of s(t; )’s are gained from the exact solution of the 
problem. It is evident that in the real-world application of this algorithm, 
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Table 2: Result of the execution of the modified RVGA, with 10 to 100 generations 
when the noise of the measurement is on the s (overspecified condition). 


Gen. | Fitness | Time(s) Ss 
10 | 1.2936 | 371.1635 | 0.1821 
20 1.1459 | 715.8062 | 0.2162 
30 | 0.6429 | 1063.0171 | 0.0823 
40 | 0.7662 | 1446.6704 | 0.0966 
50 0.4227 | 2012.2701 | 0.0515 
60 0.4109 | 2171.3517 | 0.0529 
70 0.3383 | 2543.9981 | 0.0435 
80 0.3379 | 2890.7852 | 0.0168 
90 0.2945 | 3254.4581 | 0.0080 
100 | 0.2407 | 3663.9104 | 0.0064 


these values come from sensors. It is obvious that all sensors have a small 
error of measurement. Therefore, to simulate these errors, a small random 
value is added to s(t;)’s to have noisy data. 

In this example, a(t) = 2t and 
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T(x,y,t) =e 7 * (cos(ry) + sin(ra)). 


Remark 5. In this article, we use the following equation to calculate total 
error S, in our numerical computations (see [6]): 


Ni 


l N 
= fees D(a a;)?| 
i=1 
where (N — 1) is the length of estimated vector (chromosome), @j,4 = 
1,2,---,N, come from interpolated function, and a;,i = 1,2,---,N, are 

the exact values of a(t). 


In this section, the population size of the algorithm is 20. Each chromo- 
some in the population has 1000 genes (t = 0.001, 0.002, 0.003, ..., 0.999). 
Table 2 shows the result of the execution of the modified RVGA, with 10 to 
100 generations when the noise (noisy data = input data + (0.00001)rand(1)) 
is on s. To study the ability of the algorithm to solve an IPP, we have 
severally executed it while the noise of the measurement is on the different 
functions. Table 3 illustrates the result of those executions. Also, Figures 
1-6 show the executions. 


Remark 6. One should note that by tending to the real value of a(t), the fit- 
ness value (the difference of the calculated and measured values in the sensor 
position) will be reduced (see (20)). Thus the total error S' will decrease. 
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Table 3: Result of the different executions when the noise of the measurement is on 


the different functions. 


Noisy data | Gen. | Fitness | Time(s) S) 

s(t;) and f 100 0.1381 | 3551.8039 | 0.0216 

s(t;) and go | 100 | 0.4914 | 3602.6282 | 0.0525 

s(t;) and gy 100 0.3511 | 3649.0969 | 0.0385 

s(t;) and ho | 100 0.4028 | 3704.7115 | 0.0627 

s(t;) and hy | 100 | 0.3718 | 3679.4812 | 0.0373 
2 ww i a 

: je 


Figure 1: The exact and approximated a(t) (Left) and the best chromosome (Right) 


with the noisy data on condition (22g). 
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Figure 2: The exact and approximated a(t) (Left) and the best chromosome (Right) 


with the noisy data on conditions (22g) and (22b). 


The time-dependent diffusion equation: An inverse diffusivity problem 49 


Figure 3: The exact and approximated a(t) (Left) and the best chromosome (Right) 
with the noisy data on conditions (22g) and (22c). 
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Figure 4: The exact and approximated a(t) (Left) and the best chromosome (Right) 
with the noisy data on conditions (22g) and (22d). 


8.1 Computational time and accuracy comparison 


To study computational time and accuracy efficiency, the example of this 
section is solved by a general real-valued genetic algorithm. Table 4 shows 
the execution of a general real-valued genetic algorithm for the number of 
generations from 10 to 100. Also, Figures 7, 8, and 9 illustrate a comparative 
study of execution time, fitness, and total error (5) between the modified 
and general real-valued genetic algorithm, respectively. 
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Figure 5: The exact and approximated a(t) (Left) and the best chromosome (Right) 
with the noisy data on conditions (22g) and (22e). 


“oa 05 
tis) 


Figure 6: The exact and approximated a(t) (Left) and the best chromosome (Right) 
with the noisy data on conditions (22g) and (22f). 
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Table 4: The execution of a general RVGA, with 10 to 100 generations when the noise 
of the measurement is on the s (overspecified condition). 


Gen. | Fitness | Time(s) Ss 

10 3.0323 | 318.4011 | 0.3824 
20 | 2.9706 | 672.6387 | 0.3792 
30 1.4854 | 992.8147 | 0.2140 
40 1.4103 | 1156.1270 | 0.2041 
50 1.1419 | 1434.9134 | 0.2043 
60 1.1557 | 1617.6324 | 0.1285 
70 1.0357 | 1701.0975 | 0.1253 
80 1.1049 | 1927.8542 | 0.1077 
90 0.9340 | 2154.6946 | 0.1174 
100 | 0.9267 | 2315.6700 | 0.0939 
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Figure 7: Comparative execution time study between the modified and a general 
RVGA. 
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Figure 8: Comparative fitness study between the modified and a general RVGA. 
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Total Error (S) Study 
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Figure 9: Comparative total error (S) study between the modified and a general RVGA. 


It is evident from Figure 7 that the execution time of the modified and the 
general RVGA increase steadily as the number of generations grows from 10 
to 100. Compared to the general RVGA, the modified RVGA generally has a 
higher run time due to an additional step. In contrast to the execution time, 
there are significant improvements in the fitness values and accordingly, the 
total errors (S') in the modified RVGA. This study suggests that the modified 
algorithm is more suitable for the applications that need a higher amount of 
accuracy. 


9 Discussion 


In Table 2, the amount of S (total error) significantly decreases while the 
number of generations (iterations) is increasing from 10 to 100. Our ex- 
periments showed that the amount of S fluctuates around 0.0064 after 100 
generations. Therefore, an optimal solution can be obtained in problem (1) 
with 100 generations while the noise of measurement is only on s (overspeci- 
fied condition). Table 3 illustrates the results of the algorithm when the noise 
of measurement is on one of the initial or boundary conditions in addition to 
s. It is evident that, the amount of total error increases when the number of 
noisy resources increases. Therefore, the amount of total error varies between 
0.0216 and 0.0627. 


10 Conclusion 


1. The presented modified genetic algorithm is a capable method to find 
unknown a(t) in an IPP even with noisy data. 
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2. Results showed that an acceptable estimation can be obtained by the 
presented algorithm with 100 generations within a CPU with a clock 
speed of 2.4 GHz. 


3. The presented algorithm is stable with respect to a small perturbation 
in the input data. 


4. The execution time of the presented algorithm is acceptable to be used 
in real-world applications. 
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